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We introduce and prove a novel linear response stability theory for spin glasses. The new stability 
under suitable perturbation of the equilibrium state implies the whole set of structural identities 
that characterize the spin glass phase. 



- - - The linear response stability under small perturbations is a basic method in statistical physics [1],[2]. By suitably 
^ , perturbing the energy of a system it allows to obtain relations (structural identities) among physical quantities by 
' simply computing a derivative. Those identities identify the order parameters of the theory thus reducing its apriori 
fS| , degrees of freedom. The mean field theory for the ferromagnetic phase in which the only order parameter is the 
^ magnetization is a celebrated example of such an approach. 

The purpose of this paper is to introduce a novel linear response stability for spin glasses. Their low temperature 
C/^ ' phase is still subject of intense investigations and debates in physics. The importance of the spin glass problem is not 
only related to the condensed matter original formulation but also to the numerous applications that its methods and 
^Sl techniques have brought into other disciplines [3], [4], [5]. At the core of the problem there is the much awaited solution 
[ of the debate between droplet [6] and mean field [3] pictures, as well as the trivial versus hierarchical (ultrametric) 
' I '■ organization of pure states in the low temperature phase. If such a debate is basically settled in the mean field 
O ^, approximation the physics community hasn't reached yet any consensus on the relevance of the ultrametricity for the 
' ■ short range finite dimensional models (see nevertheless the numerical results in [7]). 

] In a former linear stability theory for the spin glass phase (stochastic stability) introduced in [8] it was shown how 
' to derive an infinite family of polynomial structural identities (see [9] at equilibrium and [10] out of equilibrium). On 
^ I the other hand a larger set of identities, known as Ghirlanda-Guerra identities [11], has been obtained by the control 
of the energy fluctuations with respect to the equilibrium state first in the mean-field approximation and then also 
I for short-range models [12], [13]. 
^ ' From the linear response stability that we introduce and prove in this paper we obtain the entire set of Ghirlanda- 
00 ', Guerra identities of the spin glass theory. By consequence their derivation requires only a simple computation of 
' a derivative, expressing the stability under a suitable perturbation of the equilibrium state. The new perturbation 
reproduces both thermal and disorder fluctuations and is able to bypass the refined fiuctuation theory bounds of the 
quenched measure used by the mathematical physics community. 

The Gibbs-Boltzmann state cop^N of a statistical mechanics system of N interacting spins a ~ (cti, ctat), with 
Hamiltonian H{a) at inverse temperature /3, admits the classical probabilistic interpretation as the deformation of 
the uniform measure over spin configurations fj,N- 



OO 



>< 



^ ' with 



MA'C/) = 4E/('^)' (2) 



and / a smooth bounded function of the spin configurations. Such a deformed state ujp^N fulfills a remarkable stability 
property with respect to further small deformations (perturbations) : considering the Hamiltonian per particle 

M.) = ^ (3) 

and the perturbation with parameter A defined as 

(A) _ ^fi,N{fe-^^) 

the Gibbs-Boltzmann measure is stable, i.e. A-independent, in the thermodynamic limit N ^ oo. In fact one can 
observe that the perturbation amounts to a small temperature shift: 

(A) ... _ l^Nife-^"-^'^) _ ... ... 
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which imphcs that it has a vanishing effect in the large volume limit apart form isolated singularities, possibly related 
to phase transitions. More precisely one can prove the stability as follows: since for all /? intervals and all values of A 
one has, thanks to (5), 



one obtains: 



4uNif)-4'J,Nif) 

N 



ffJi duj'-^^ (f) 

hm / ^f' dp = V A, V [/3o,/3i] . (8) 

As a consequence, computing the derivative at A = 0, one has 

lim / Kiv(//i)-w^,Jv(/V/3,iv(ft)]d/3 = 0. (9) 

For the special case f = h the previous formula implies that the Hamiltonian per particle converges to a constant for 
large volumes with respect to the Gibbs measure, at least in /3— integral average. Higher order derivatives with respect 
to A of (i.e. cumulants of h) are then enforced to vanish since cumulants are homogeneous polynomials of 

the constant values of h with coefficients whose sum is zero. 

Formula (9) has interesting consequences. It says, for instance, that the order parameter (i.e. the magnetization) 
for a mean field ferromagnetic Hamiltonian has a trivial distribution. In the Curie- Weiss model at zero magnetic field, 
for which the Hamiltonian per particle is the square magnetization, the previous identity implies (by choosing f = h) 
that in the thermodynamic limit, 

W/3(o-i (72 0-30-4) = a;/3(0-lO2)^ (10) 

in /3— average [14]. One can indeed prove that (10) holds for all ^ using the methods developed in [15]. The choice 
/ = /i", or cquivalcntly higher order derivatives in A of the perturbed state, gives the well known factorization property 
of the 2n-point function as an n-th power of the 2-point function. 

In a disordered system defined by a centered Gaussian Hamiltonian y,{a) of covariance (generalized overlap) 

AvCH(cr)-H(r)) ^ Ncn{(t,t) (11) 

the equilibrium measure is the quenched average of the random Boltzmann-Gibbs state LUp^M'- for a bounded random 
function / it is defined by 

{f)p,N = kY{LOp,N{f)) . (12) 

The thermodynamic properties of the system are expressed in terms of a set of random variables {cjj } related 
to the quenched expectation of the covariance entries. Namely, considering the Boltzmann-Gibbs product state 
^fi,N = i^fi,N X i^fi,N, one defines the random variables Cij and their joint distribution by: 



dl3 

(7) 



Ef. 



idj) = Av (0^,;v(c(aW , a(^)))) . (13) 



In [8] it was identified a stochastic stability property of the quenched state, i.e. an invariance with respect to the 
stochastic perturbation: 

where /C(o) is a random field (independent from the Hamiltonian) whose covariance is CAr(o, r), and it was shown 
that the stochastic perturbation is equivalent to a temperature shift 
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from which stabiUty follows [12] 



lim / = V A, V[/3o,/3i]. (16) 



By consequence at A = one obtains: 

lim / Av(u;^,jv(A)-W/3,Jv(/)w^,iv(/i))rf/3 = 0. (17) 



AT 

The previous formula implies (taking f — h and integrating by parts) 

im / 

N 



lim / E;3,;v(Ci2 -4C1,2C2,3 + 3C1,2C3,4) = . (18) 

We stress the fact that the previous identity holds for a general Gaussian Hamiltonian, both mean field or short range, 
in terms of its own covariance. For / = /i" one can see [12] that the identities that can be derived from (17) are, like 
the (18), zero average polynomials in the Cij with respect to the quenched measure. See also [9] for an alternative 
derivation. 

In [16] it was introduced a method, based on bounds for the energy fluctuations, which leads to the set of Ghirlanda- 
Guerra identities [11, 13, 17, 18]; the lowest order are for instance: 

E/3,JV(C1,2C2,3) = ^E/3,W(C?,2) + ^%,W(C1,2)' (19) 



%,JV(C1,2C3,4) = ^ %,JV(C?,2) + I %,JV(C1,2)' • (20) 

Unlike the set of identities that can be derived from stochastic stability, these also include non-linear terms of the 
overlap expectations. In recent times it was shown that an invariance under reshuffling introduced in the framework 
of competing particle systems [19] implies the whole set of Ghirlanda-Guerra identities [20]. 

To this purpose we introduce and prove here a novel stability property for the spin glass quenched state. We show 
that from such a stability property the whole set of Ghirlanda-Guerra relations can be derived. 

We define the perturbation of quenched state as 



//.u(A) _ Av (a.^,^(/e-^'')) 



We observe that this new perturbation is the analog, for the quenched measure of a random Hamiltonian, of the 
standard perturbation (4) introduced for deterministic systems with respect to the Boltzmann-Gibbs measure. On 
the other hand we notice that while the stochastic stability perturbation (14), as much as the standard perturbation 
for deterministic system, amounts to a small temperature shift, the newly introduced perturbation cannot be reduced 
to just a small temperature change (see the remark after formula (32)). More precisely the explicit expression of (21) 
reads 



Av I "j^^ e-fHi") 



where it clearly appears that only the numerator of the random Boltzmann-Gibbs state is affected by the change. 
Our main result is summarized by the following 

Proposition With the definition given above, the quenched state of a Gaussian spin glass is stable under the defor- 
mation (21), i.e. 



lim / 



hm ' --'^'^ 



N->-00 I a dX 



d/3 = Q. (23) 

A=0 
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Moreover the property (23) implies the whole set of the Ghirlanda-Guerra identities: for a hounded f function of the 
generalized overlaps {cjj} (with i,j € {1, ■■■,n}) the quantity 5Ar(/3) defined by 



lE/3,Jv(/ci,„+i) = -E;3,jv(/)]E;3,jv(ci,2) 



1 " 

-^E^,^(/ci,,)+(5Ar(/3) 



(24) 



j=2 



goes to zero in ^-average and for large volumes. 
Proof: A simple calculation shows that 



dX 



= {fh)^,N- {f)0,N{h)^,N . 



(25) 



A=0 



The right hand side can be decomposed into two terms which can be identified as the thermal and the disorder 
correlations: 



d{{f))f^N 



Av (uj^^Nifh) - UJi3^N{f)0J^,N{h)) + 



A=0 



+Av (a;^,iv(/)w,3,jv(/i)) - Av {iJ0,N{f)) Av {w^,N{h)) . 



(26) 



In [13] the two previous terms were proved to converge to zero in /3 average and using integration Gaussian by parts 
it was shown how they imply formula (24). □ 

Remarks 

i) It is interesting to notice that the new stability property introduced in this paper as well as those introduced in the 
past admit a simple formulation in terms of cumulant generating function. Defining that function for the quenched 
state as 



^;.,;v(A) =lnAv(^^) = ln(e^^), 
the (23) is equivalent, by the boundedness of f and the Schwarz inequality 



\{fh) {f){h)\ < ^)^{h^)-{h}^ 

to the property of asymptotic flatness at the origin 



lim 

N->-fx 



rPi 
ho 



dA2 



d/3 = 0. 



A=0 



In particular defining the generating function of thermal fluctuations as 

and the generating function of disorder fluctuations as 

^^_jv(A) = lnAv(e^"''-"('')) 

one has 



dX^ 



(f'4)l3,N{X) 



A=0 



rfA2 



+ 



A=0 



rfA2 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



A=0 



which shows that the new perturbation can he seen as a sum of thermal and disorder perturbations. 
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ii) The statement of the Proposition can be made to hold point-wise in /3 under suitable technical assumptions, 
namely the differentiability of the free energy. This indeed has been shown for the standard stochastic stability in [21 J 
and for the Ghirlanda-Guerra identities in [22]. 

The results shown in this paper provides a straightforward method to obtain the structural identities of the spin 
glass phase known as Ghirlanda-Guerra identities by a simple computation of a derivative and a Gaussian integration 
by parts. This provides a new interpretation, using a stability argument, of the vanishing fluctuation property from 
which they were originally derived [11]. The relevance of the stability properties and of the Ghirlanda-Guerra identities 
has been shown in the work [19] and [23] where, under the hypothesis of discreteness of the overlap distribution it was 
proved. rcspc;ctivc;Iy. that competing particle systems satisfying invariance under reshuffling or spin systems satisfying 
Ghirlanda-Guerra identities do fulfill the hierarchical structure (ultrametricity) originally introduced in the Parisi work 
for the mean field spin glass [3]. The present work provides a further bridge between those two approaches, whose 
mutual relation has still to be fully clarified [24] , suggesting that the invariance under reshuffling is well represented 
by our newly introduced stability under perturbation. 
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